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ABSTRACT 

The formation and growth of galaxy disks over cosmic time is crucial to our understanding of galaxy 
formation. Despite steady improvements in the size and quality of disk samples over the last decade, 
many aspects of galaxy disk evolution remain unclear. Using two square degrees of deep, wide-field 
i'-band imaging from the Canada-France-Hawaii Telescope Legacy Survey, we compute size functions 
for 6000 disks from z = 0.2toz = l and explore luminosity and number density evolution scenarios 
with an emphasis on the importance of selection effects on the interpretation of the data. We also 
compute the size function of a very large sample of disks from the Sloan Digital Sky Survey to use as 
a local {z ~ 0.1) comparison. CFHTLS size functions computed with the same fixed luminosity-size 
selection window at all redshifts exhibit evolution that appears to be best modelled by a pure number 
density evolution. The z = 0.3 size function is an excellent match to the z = 0.9 one if disks at the 
highest redshift are a factor of 2.5 more abundant than in the local universe. The SDSS size function 
would also match the z = 0.9 CFHTLS size function very well with a similar change in number 
density. On the other hand, the CFHTLS size functions computed with a varying luminosity-size 
selection window with redshift remain constant if the selection window is shifted by 1.0—1.5 mag 
towards fainter magnitudes with decreasing redshift. There is a weak dependence on disk scale length 
with smaller {h < A kpc) disks requiring more luminosity evolution than larger ones. Given that 
changes in number density are primarily due to mergers and that current estimates of merger rates 
below z = 1 are low, luminosity evolution appears to be a more plausible scenario to explain the 
observations. 

Subject headings: galaxies: formation, evolution 
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1. INTRODUCTION 

The size evolution of disk galaxies provides a unique 
constraint for galaxy formation models. In the ACDM 
paradigm, structure forms hierarchically. Primordial 
density perturbations gradually attract dark matter 
into haloes that become the sites of galaxy forma- 
tion. Through tidal torques, these haloes and their gas 
contents acquire angular momentum. The gas within 
the halo cools and condenses to form a flattened disk, 
supporting itself through rotation (Fall and Efstathiou 
1980). The first disks formed in these haloes are small 
and dense. After several mergers as well as gas accretion, 
these disks grow to become giant spirals and subsquently, 
ellipticals (Steinmetz and Navarro 2002; Toomre and 
Toomre 1972). Mergers are thus an important mecha- 
nism in the evolution of galaxies. They are also quite 
violent processes. Minor mergers are thought to lead to 
the growth of a bulge (Barnes and Hernquist 1992) and 
to an increase in disk thickness through heating (Toth 
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and Ostriker 1992). Equal mass mergers can completely 
change the original morphology of a galaxy by turning 
disks into spheroids. Although the gas content of the 
constituent galaxies may permit a re-growth of a disk 
around a bulge galaxy (Robertson et al. 2006), a pure, 
bulgeless disk galaxy grows mainly through quiescent ac- 
cretion. The cold, fragile structure of disk galaxies thus 
make them unique probes of hierarchical mass assembly; 
they provide a benchmark for the sizes of galaxies at any 
epoch as they give the sizes of objects that have been 
relatively undisturbed through cosmic time. 

Generally, models predict that large disk galaxies 
should only be stable at recent {z < 1) epochs (Mo et al. 
1998; van den Bosch 1998). However, previous surveys 
(e.g., Lilly et al. 1998; Roche et al. 1998) have found 
that the size distribution of disks remained relatively un- 
changed at z < 1 once about 1 mag of luminosity evo- 
lution had been taken into account. This suggests that 
size evolution and the bulk of mass assembly must have 
occurred prior to z — 1. Although luminosity evolu- 
tion is a natural consequence of an aging stellar popu- 
lation (Tinsley 1978), previous surveys have also found 
a range of values; from little or no evolution (Simard 
et al. 1999; Ravindranath et al. 2004) to as much as 1 
magnitude or more (Schade et al. 1996; Lilly et al. 1998; 
Harden et al. 2005; Sargent et al. 2007). One critical as- 
pect of these surveys is the treatment of selection effects 
as different treatments can lead to different conclusions 
from the same dataset. Despite concluding otherwise, 
both Simard et al. (1999) and Ravindranath et al. (2004) 
were able to produce a 1 magnitude evolution with their 
respective samples by redefining their selection criteria. 
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This is a systematic problem that is not alleviated with 
larger samples. For example, one may use a fixed sur- 
face brightness survey window with redshift and follow 
changes in galaxy properties within this window. Al- 
ternatively, one may attempt to track the same galaxy 
population across all redshifts by sliding the survey win- 
dow according to some assumptions on the evolution of 
this population. 

Earlier disk galaxy surveys did not have the advan- 
tage of the wide imaging capabilities afforded by recent 
studies. These earlier observations had to deal with the 
additional complexity of cosmic variance due to sam- 
pling in small fields. As survey field sizes have grown 
larger, sample sizes have also grown substantially. As 
a result, more statistically significant samples are now 
available, and the expanded sky coverage of these sur- 
veys minimizes the effect of cosmic variance. Over a 
decade ago, Vogt et al. (1996) used a sample of 16 galax- 
ies and found that large, massive disks were in place by 
z K 1. Roche et al. (1998) used a more sizeable sample 
of 347 objects and found that disks have undergone a 
surface brightness evolution of approximately 0.95 mag 
over 0.2 < z < 0.9 through a combination of size and 
luminosity evolution with most of the size evolution oc- 
curring prior to = 1. Using 341 objects selected from 
the Canada-Prance-Hawaii Redshift survey, Lilly et al. 
(1998) fomid that the size function of larger disks (disk 
scale length h > 'S.2h^^ kiloparsecs) is approximately 
constant to z w 1 and that disks are approximately 0.8 
mag brighter at z = 0.7. These studies seem to sup- 
port a model where large disk galaxies are in place by 
z = 1 in the same number density and physical sizes as 
they are today. In addition, they claim that the luminos- 
ity evolution is consistent with a simple, passively aging 
population. However, other studies find that the lumi- 
nosity evolution in disks is minimal over this time. With 
a sample of 190 field galaxies from the Deep Extragalactic 
Evolutionary Probe (DEEP), Simard et al. (1999) per- 
formed a careful analysis of selection effects and found 
that there has been no evolution in disk surface bright- 
ness. Ravindranath et al. (2004) used 1508 galaxies from 
the Great Observatories Origins Deep Survey (GOODS) 
and also found that there has been little or no evolution 
(< 0.4 mag) in the surface brightness of disks over the 
redshift range 0.2 < z < 1.25. More recent studies came 
to the conclusion that disk-dominated galaxies actually 
evolved in luminosity. Barden et al. (2005) used 5664 
disk-dominated galaxies from the Galaxy Evolution from 
Morphologies and SEDs (GEMS) survey and found that 
there has been a brightening of approximately 1 magni- 
tude to a; < 1.1. By applying the same amount of surface 
brightness evolution as measured by Barden et al. (2005) 
to a sample of 36 disks in the Hubble Ultra Deep Field, 
Trujillo and Pohlen (2005) find that these disks have 
had moderate (~ 25%) inside-out growth since z = 1.1. 
As part of the COSMOS project, Sargent et al. (2007) 
examined approximately 12,000 disk-dominated galaxies 
and found that the number density of large disks has re- 
mained fairly constant out to z = 1. Of these however, 
they found that their very largest disks {h > 10 kpc) are 
only 60% as abundant as today, but inferred that some 
of these objects may have a more significant bulge com- 
ponent in the present day. This could be interpreted as 
a transformation from bulgeless disks at z = 1 into rel- 



atives of the Milky Way. Studies at high redshift now 
benefit from a better imderstanding of local galaxy pop- 
ulations thanks to the Sloan Digital Sky Survey (Strauss 
et al. 2002), but a clear consensus on the evolution of 
disk galaxies since z = 1 has yet to emerge. This pa- 
per focuses on the size evolution of galaxy disks out to 
~ 1 from the deep, wide-field imaging of the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS). It is 
organized as follows. The data and bulge-|-disk decom- 
positions are described in Sections 2 and 3. The sample 
selection and selection effects are discussed in Section 4. 
In particular, we will explore how different treatments 
of selection effects may be interpreted as different evolu- 
tionary models. The disk size functions are presented 
in Sections 5 and 6. The implications of our results 
for the evolution of galaxy disks are discussed in Sec- 
tion 7. The cosmology adopted throughout this paper is 
{Ho,n^,nA) = (70, 0.3 , 0.7). 

2. IMAGING DATA 

The images used for this study were obtained as part of 
the Deep component of the Canada- France-Hawaii Tele- 
scope Legacy Survey (CFHTLS-Deep). The CFHTLS is 
a large five-year project that began in 2003. The overall 
goal of CFHTLS-Deep is to obtain a better understand- 
ing of the early universe through detection of thousands 
of supernovae and the study of the galaxy population. 
Covering nearly four square degrees, the survey uses 
large statistical samples to build stronger constraints on 
galaxy evolution and the star formation history of the 
Universe. Such a large survey area is accessible thanks 
to the MegaPrime prime-focus mosaic imager (9x4 mo- 
saic of 2080 pixel x 4622 pixel detectors). With a scale 
of 0.187" /pixel, the full mosaic has a 0.96 deg x 0.96 
deg field of view. The morphological analysis was per- 
formed using the i' band images, and other bandpasses 
were used for the computation of photometric redshifts. 
We used two deep fields for this study: Dl (a= 02:26:00, 
S =-04:30:00) and D3 (q:=14:17:54,+52:30:31). The D3 
field partially overlaps the well-known Groth-Westphal 
Strip region (Rhodes et al. 2000). Analysis for each field 
was performed on stacked i'-band images (6580-9120A) 
consisting of 150 exposures totaling 12.9 hours and 96 
exposures totaling 13.9 hours, respectively. The depths 
of our Dl and D3 image stacks (5-cr detection limit) were 
26.5 ± 0.1 mag and 26.7 ± 0.1 mag respectively. 

All previous studies of disks at high redshift have used 
space-based imaging from the Hubble Space Telescope. 
While HST imaging clearly offers better spatial resolu- 
tions, ground-based observations with wide-field imagers 
can deliver images over larger fields with longer integra- 
tion times. Large field size is important to reduce the 
effects of cosmic variance which can be significant even 
at square degree scales (Somerville et al. 2004; Trenti and 
Stiavelli 2007). Longer integration times mean that lower 
surface brightness levels can be probed. The limited spa- 
tial resolution of ground-based observations is offset by 
an increase in the radial range over which photometric 
decompositions can be performed. While it is true that 
the structural parameters of high-redshift bulges cannot 
be measured from the ground due to seeing-limited reso- 
lution, disk scale length measurements can compare very 
well with space-based results thanks to larger disk ap- 
parent sizes and larger fitting ranges. 
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3. METHODOLOGY 
3.1. Image Processing and Point-Spread-Function 

Basic preprocessing (bias removal, fiat fielding and 
fringe correction) was performed using the Elixir pipeline 
(Magnier & Cuillandre 2004), and elixir-processed im- 
ages were then retrieved from the Canadian Astronom- 
ical Data Centre. All the images were astromotrically 
and photometrically calibrated and then stacked nsing an 
early version of the MegaPipe reduction pipeline (Gwyn 
2008). Sources are extracted using the SExtractor rou- 
tine (Bertin and Arnouts 1996), and a catalogue of ob- 
jects is generated. Basic photometry in all 5 bands was 
done using Kron-style apertures using the double im- 
age mode of SExtractor. Photometric redshifts and 
rest-frame absolute magnitudes were computed follow- 
ing Gwyn (2008). Photometric redshifts were calculated 
using a template fitting method. In this method, galaxy 
broadband colours in Megacam u*, g', r', i' and z' are 
used to create a low resolution spectral energy distri- 
bution (SED) which is then compared to a number of 
known galaxy template spectra (eg. star-forming, quies- 
cent etc). Here we use template spectra from Coleman 
et al. (1980). For a single galaxy, a suitable template is 
determined by the one that returns the lowest value. 
Using features that are common to both spectra, the 
galaxy spectrum is shifted until it is best aligned with 
the template distribution. This shift is used in the stan- 
dard redshift relation z = (A — Ao)/Ao to determine the 
redshift. There is not a single "feature" that is used for 
Ao; the optimal shift is found by an overall agreement 
between the galaxy SED and the template. The relative 
error in the photometric redshift measurements is about 
Sz/{1 + z)^ 0.11 (Nuijten ct al. 2005). 

The point spread function of the images was built using 
the DAOPHOT (Stetson 1987) package within the Im- 
age Reduction and Analysis Facility (IRAF). A variable 
Moffat point-spread function (PSF) was constructed for 
each chip. Typically, 75 point source objects were used to 
generate the model point spread function for each chip. 
Selected objects were in focus, had very high SExtractor 
stellarity indexes and were fairly well isolated with regu- 
lar, circular isophotes. A model for the point spread func- 
tion was determined through several iterations. Point 
sources were convolved with the PSF model and then 
subtracted from the original image at the locations of 
the objects used to derive the PSF. The subtracted im- 
age was then visually inspected for residuals. A good 
PSF model will yield clean image residuals with no sys- 
tematics present. Possible systematic errors for the point 
spread function include a very bright or dark spot in the 
center of the object or dark spots within the PSF fit- 
ting radius. Parameters were then adjusted based on the 
type of systematics present, and the model was recom- 
puted. Once all systematics were minimized, the point 
spread function was tested with our surface brightness 
fitting routine (see section 3.2) and minor adjustments 
were made when necessary. 

The typical full-width-half maximum of the stacked 
images was 0'.'95, although it had to be adjusted from 
chip-to-chip. Chips further from the center of the mo- 
saic generally had larger FWHM values than those near 
the center, but they never deviated by more than 10% 
from this value. On a given chip, the variability of the 



point spread function was also modelled by DAOPHOT. 
A 'true' point spread function was derived for the entire 
chip, but several lookup tables were computed to mea- 
sure deviations from this model. The point spread func- 
tion was allowed to vary quadratrically over the chip and 
thus included terms proportional to l,x, y, , xy and y^ 
where x and y are positions on the chip. The choice 
of the Moffat function to fit the stellar profile is fairly 
common. It is known to be numerically well behaved in 
fitting narrow point spread functions, and its wings also 
fit stellar profiles better than a gaussian profile (Trujillo 
and Agucrri 2004). The PSF fitting radius was 5'.'6-6'.'0. 
Larger values were found to contaminate the point spread 
function with neighbours, and smaller values did not pro- 
vide a good model of the innermost pixels. 

3.2. Bulge + Disk Decompositions 

Measuring galaxy structural parameter from 
bulge-|-disk (B-|-D) image decompositions has been 
used extensively in the past (Simien and de Vaucouleurs 
1986; Kent 1985) and continues into the present day 
(Marlcau and Simard 1998; Peng ct al. 2002a). Here, 
we use the surface brightness fitting routine "Galfit" 
described in detail in Schade et al. (1995) and Saintonge 
et al. (2005). This is not to be confused with another 
routine developed for structural decomposition of galaxy 
images also called 'GALFIT' and developed by Peng 
et al. (2002b). For each of the 400,000 galaxies in 
our sample. ''Galfit'' creates a 'postage stamp' of the 
object from the stacked image based on its SExtractor 
coordinates. Upon finding an accurate centroid for the 
object, a 'symmetrized' image is generated in order 
to eliminate light contamination from neighbouring 
objects. Symmetrization is achieved by taking the 
original postage stamp, rotating it by 180°, subtracting 
this rotated frame from the original and only preserving 
features that are at least 2(7 above the noise level of 
the image. Only these features that remain after the 
symmetrization process are used in the fitting process. 
The sky background level is not fit; an estimate of the 
background is input into the routine and then subtracted 
from the postage^ stamp. The symmetrized object used 
by Galfit is seen in Figures 1 and 2 in the top row of 
images in the second panel from the left hand side. A 
detailed analysis of possible photometric biases due to 
symmetrization is given in Saintonge et al. (2005). 

The routine convolves the point spread function with 
three idealized surface brightness models (pure r^/^ 
bulge, pure disk and bulge-|-disk) and computes a 
for each. The routine fits 6 or 7 parameters for each 
of the idealized models. The actual set /number of fit- 
ting parameters depends on the particular model. For 
example, parameters for the bulge+disk model are bulge 
effective radius, bulge axial ratio, bulge position angle, 
disk scale length, disk axial ratio, disk position angle and 
fraction of bulge light to total galaxy light {B/T). The 
routine begins with preset initial parameter values and 
varies each parameter in the direction which decreases 
the of the fit. Initial parameters are globally set to 
some intermediate value (e.g., B/T always begins at 0.5, 
position angles at 45 degrees). The routine is sensitive 
to such input values only for objects that are below the 
detection limit of the survey; otherwise it generally con- 
verges to the same best-fit model regardless of initial val- 
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ues. The global minimum is determined by a modified 
Levenberg-Marquardt algorithm. 




Fig. 1. — Example of a bulgc+disk fit. Top row from left to right: 
Postage stamp galaxy image extracted directly from science image, 
symmetrized galaxy image with neighbours removed, and residual 
images from the bulge, disk and bulgc+disk model fits. Middle row: 
The values for the pure bulge, pure disk and bulge+disk models 
are 5.00, 1.38 and 1.21 respectively. Bottom row. Bulge and disk 
components of the bulge+disk model. In this case, as is true for 
the disks in our sample, the bulge component of the bulge+disk 
model is a small fraction of the total light for the galaxy. This 
galaxy has a bulge fraction of 9%. 




Fig. 2. — Example of a pure disk fit. The layout of the images is 
same as Figure 1. This galaxy was best fit by a pure disk model. 
As this image illustrates, the bulge+disk model (x^=1.88) con- 
verged to a pure disk model (x^=1.86). (Bulge model: X'^= 6.56) 
Although in this case the disk model was selected, both models 
yielded the same parameters for the galaxy, and the values be- 
tween the models differed by less than 1%. 



4. SAMPLE SELECTION 

4.1. Disk selection 

For every galaxy in the sample, Galfit returns val- 
ues for the pure bulge, pure disk and bulge+disk model 
fits to the image. Objects that did not converge to a 
solution were not used in the sample described below. 
Approximately 4.1% of the initial set of objects did not 
converge. Many of these were objects whose neighbours 
could not be properly removed or had other structural 
peculiarities. Some of these are discussed in Saintonge 
et al. (2005). We selected galaxies that had: 1) the low- 
est produced by a pure disk model or 2) the lowest x^ 
produced by a bulge-|-disk model and B/T < 0.2. The re- 
sulting sample includes nearly 65,000 galaxies plotted by 
redshift in Figure 3. Our disk selection must also include 
a lower cut on disk size to exclude galaxies with unreli- 
able disk scale lengths that may affect the size function. 
In order to determine the smallest apparent galaxy size 



that could be reliably measured, we used a set of 5000 
artificially generated galaxy images. 

The input objects were bulge-|-disk models with vary- 
ing bulge fractions. The artificial galaxies were convolved 
with the real point spread function and then inserted 
into real image stacks. They were then analyzed with 
Galfit in exactly the same way as real galaxies with pure 
bulge, pure disk and bulge-|-disk models. We first com- 
pared the number of objects returned by Galfit versus 
the number of input objects as a function of the input 
disk scale length. We found that the number of simu- 
lated galaxies with recovered disk scale lengths less than 
0'.'37 was as high as 10 times in excess of the input num- 
ber. We carefully examined many "problem" objects to 
visually confirm the effect. The effect was particularly 
problematic for sizes smaller than approximately 0'.'21. 
For sizes between 0'.'21 and 0'.'32, the fraction of objects 
returned improved drastically; Galfit would return a very 
reasonable 1.1-1.5 times the number of input objects at 
0'.'32. However, at sizes greater than 0'.'37 the number re- 
turned matched the input number exactly. We therefore 
restricted our final disk sample to disks larger than this 
threshold to obtain a reliable size function. Our mini- 
mum threshold disk scale length corresponds to 3 kpc at 
z = \. We also used simulations to further characterize 
the region in size and luminosity over which disks could 
be reliably detected and measured out to z = 1. Given 
our magnitude limit of i'=24.5, surface brightness dim- 
ming with increasing redshift pushes galaxies out of our 
sample, and larger disks are preferentially excluded be- 
cause their light is spread out over a greater area. Our 
luminosity and size selection criteria are shown in Figure 
3. 

4.2. Selection Effects 

In order to evaluate the importance of selection effects 
on the calculation of the disk size function, the bound- 
aries of the survey selection window must be carefully 
characterized. The treatment of selection effects can be 
done by adopting two different approaches. These ap- 
proaches are not at all equivalent- though they have been 
presented as such in previous studies of high-redshift 
disks. One approach is to determine the selection window 
of a survey at its high redshift end (where selection effects 
are expected to be most severe) and to apply this fixed 
selection window at all redshifts to map changes in galaxy 
properties within it. The other approach is to attempt to 
track the evolution of a given galaxy population across 
all redshifts by using a selection window that is varying 
with redshift based on a set of assumptions (e.g., pure 
luminosity evolution) regarding the expected evolution 
of the population of interest. An illustrative analogy to 
these two approaches is the two ways in which one can 
formulate the equations of fluid dynamics. The "Eule- 
rian" formulation uses coordinates that describe what 
happens at a fixed point in a fluid volume whereas the 
"Lagrangian" formulation uses coordinates that follows 
a given fluid volume element as it moves around space. 
The "Eulerian" (e.g., Simard et al. 1999; Ravindranath 
et al. 2004) and "Lagrangian" (e.g., Schade et al. 1996; 
Harden et al. 2005; Sargent et al. 2007) treatments of se- 
lection effects in galaxy surveys have both been used in 
previous studies even though they have been presented 
using the same language, and this ambiguity may be re- 
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sponsible for previous studies yielding results that appear 
discrepant at first sight. 

Here we will use both approaches to the treatment of 
selection effects in our survey to illustrate how important 
the choice of approach is to the interpretation of the re- 
sults from the same set of data. Unfortunately, neither 
approach is perfect. The "Eulerian" approach may be 
tracking different populations as they move in and out 
of the fixed visibility window whereas the "Lagrangian" 
approach may mistake trends in a non-evolving popu- 
lation for real evolution. This last point is particularly 
important. For example, local galaxy disks do not have 
a constant central surface brightness (de Jong and Lacey 
2000; Shen et al. 2003; Driver et al. 2005, Simard 2008, 
in prep.), and the surface brightness of bright disks (the 
ones visible at z = 1) can be up to 1 magnitude brighter 
than the surface brightness of fainter disks. Such a vari- 
ation in surface brightness could be interpreted as evolu- 
tion by observing the same population at different red- 
shifts with a selection window that does not look at the 
same range of disk luminosities. 

To treat selection effects with a fixed visibility win- 
dow, we begin by constraining the sample to the lim- 
its defined in the highest redshift bin. These limits are 
then applied to all lower z bins as illustrated in Figure 
3. Galaxies in the highest redshift bin (0.8 < z < 1.0) 
are the most constrained in luminosity-size space with 
respect to the other redshift bins as only the brightest 
galaxies remain detectable. Imposing the constraints in 
luminosity and size as a function of luminosity (surface 
brightness) present at high redshift on every redshift bin 
will ensure that the sample used here is free of surface 
brightness selection biases. This constraint ensures that 
comparable types of galaxies are used across all redshifts 
considered. The size function is computed only with the 
galaxies in this window. 

Applying the same 'window' from the highest z bin to 
all lower bins assumes no surface brightness evolution in 
the computation of the size function. In order to evaluate 
the effects of luminosity evolution on the size function, 
this window was moved by different amounts in mag- 
nitude as a function of redshift, and the size function 
was recomputed in each redshift bin. Different amounts 
of evolution were explored, and size functions were re- 
computed as a function of luminosity evolution. Various 
surface brightness models are plotted in Figure 3. They 
will be discussed later in section 6.2. 

5. THE SIZE FUNCTION 

The size function (Schade et al. 1995; Lilly et al. 1998) 
is a powerful tool to characterize the size evolution of 
galaxies. The size function gives the space density of 
disks of a given scale length at a given redshift. It is 
important to emphasize that all previous works with the 
exception of Lilly et al. (1998) and Sargent et al. (2007) 
have focussed on the zeropoint of the luminosity-size re- 
lation of high-rcdshift disks. Therefore, they could not 
detect any evolution in luminosity and/or size that might 
have resulted in the same zeropoint. The size function 
provides a more stringent measure of disk evolution be- 
cause it also places constraints on the density (i.e., dis- 
tribution) of disks in luminosity and size. 

Galaxies entering the calculation of the size function 
must first be weighted according to how well they are 
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Fig. 3. — Selected objects in disk luminosity-size space. Dot- 
dashed red lines: selection window used to compute the size func- 
tion for the no-evolution case. Dashed coloured lines: the 'shifted' 
windows used to compute the size function with luminosity evolu- 
tion. In the 0.6 < z < 0.8 and 0.4 < z < 0.6 bins, the green line 
indicates 0.2 mag of evolution, and the blue line indicates 0.4 mag. 
The 0.2 < 2 < 0.4 bin shows the cases of 0.5, 1.0 and 1.5 mag of 
evolution shown in green, blue and violet, respectively. 

represented in the sample on the basis of their luminos- 
ity and redshift. Bright galaxies will be visible over larger 
cosmological volumes and would overwhelm galaxy sam- 
ples if they were not also less abundant than faint ones. 
The appropriate weighted number densities are calcu- 
lated using the l/Vmax method as described in Schmidt 
(1968). The accessible (or maximum) volume of a galaxy 
is defined as the comoving cosmological volume that a 
galaxy may occupy while remaining within the selection 
criteria of the survey. A galaxy of a given luminosity 
can be placed anywhere within this volume and still be 
detected. For the sake of simplicity, Vmax was computed 
with luminosity as the limiting factor as described below. 
Although the size of the galaxy may also determine the 
value of 2;ma2;(Shen et al. 2003), we find here that lumi- 
nosity is the limiting factor in greater than 99% of cases 
and thus did not consider size limitations to Vmax- To 
compute Vmax , the sample is subdivided into a number of 
redshift shells. Let Zmin,s and Zmax,s be the redshift lim- 
its of a shell, and let Zmin,m and Zmax,m be the redshift 
limits between which the galaxy meets the magnitude se- 
lection of the sample. Then the lower and upper redshift 
bounds for computing the maximum volume of a galaxy 

are Zmin — niax(z^j^ g , Zmin,m,^ and Zmax — niin(z^Qa; s 

, Zmax,m.)- One Can integrate the comoving volume ele- 
ment between these redshift limits for each galaxy: 



Vmax = ^ / dnf{0,i 



DUz) 



where Dh is the Hubble distance c/Hq, Hq is the Hub- 
ble constant, flM is the parameter for matter density in 
the universe, and Oa is the cosmological constant (Hogg 
1999). The survey area Qs corresponding to our two 
CFHTLS-Deep fields is 1.83 square degrees. The sam- 
pling fraction f{9, (f>) as a function of position on the sky 
is constant over the survey area. Dc{z) is the comoving 
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Fig. 4. — The size function for disks with h > 3 kpc for the no 
evolution case. Galaxies included here reside within the 0.8 < z < 
1.0 selection window of Figure 3. The 0.8 < z < 1.0 size function 
is shown as a red line in all the redshift bins for comparison. 

distance defined as 



Dc{z)^Dh 



dz' 



(2) 



Once 1/Vmax corrections are available for all galaxies 
in the sample, the space density of galaxies with disk 
scale lengths in the range h,h + dh at redshift z is then 
given by the sum: 



N 1 

^{h,z)dh^Y.T^ 



(3) 



where h is the scale length of the disk, and N is the num- 
ber of disks in the sample at z. The size bin widths used 
here are 0.5 kiloparsecs. and size functions are computed 
in Az = 0.2 redshift shells ( 0.2 < z < 0.4, 0.4 < z < 0.6, 
0.6 < z < 0.8 and 0.8 < z < 1.0). All error bars shown 
in the size functions presented here are computed us- 
ing bootstrap resampling. The size function is computed 
1000 times, each time using a different sample drawn 
from the original sample. The error bars represent 99% 
confidence intervals based on the distribution of values 
from all the size function realizations, and they show the 
stability of the size function values within the data set. 

6. RESULTS 

6.1. Size Function with Fixed Selection Window 

The size function of disks as a function of redshift for 
the no-evolution case is shown in Figure 4. It was com- 
puted by applying the same selection window at all red- 
shift. The selection window is the one shown for the 
redshift range 0.8 < z < 1.0 in Figure 3. As mentioned 
earlier, this window places the most stringent limitations 
on size and luminosity of the galaxies used to calculate 
the size function. In this case, the size function is not 
tracking changes in a given disk population. Rather, it 
is measuring changes in the sizes of the disks that popu- 
late the fixed selection window at different redshifts. 

The shape of the size function in Figure 4 is quite sim- 
ilar in all redshift slices. This suggests that the relative 



distribution of sizes is essentially constant over this red- 
shift range i.e., there is no size-dependent size evolution. 
The size functions of the two highest redshift bins seem 
to be in good agreement, but there is a significant differ- 
ence between the size function in the lowest and highest 
redshift bins. This difference appears as a change in the 
normalization of the size function. Such a change may 
at first glance be indicative of a change in the number 
density of disks with redshift; however Figure 5 shows 
that a pure number density evolution model would re- 
quire about 2.5 times more disks at 0.8 < z < 1.0 than 
at 0.2 < z < 0.4 and 1.5 times more disk at 0.8 < z < 1.0 
than at 0.4 < z < 0.6 in order to reproduce the observed 
size functions. The plausibility of such a model is dis- 
cussed in Section 7.2. 
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Fig. 5. — Size function of disk galaxies with number density 
evolution in the lowest two redshift bins. The Size function at 
0.2 < z < 0.4 with a multiplicative factor of 2.5 applied is shown 
in black. Size function at 0.4 < z < 0.6 with a factor of 1.5 applied 
is also shown in black. The red line shows the size function at 
0.8 < z < 1.0 bin and is reproduced in all other redshift bins. 
Points in grey show the no-evolution size functions from Figure 4. 



6.2. Size Function with Varying Selection Window 

In order to compute the size function with a varying 
selection window, we first specify our constraints in size- 
magnitude space at 0.8 < z < 1.0 as described earlier, 
and subsequently slide this window towards fainter ab- 
solute magnitudes with decreasing redshift until the size 
function at each redshift matches the high redshift one. 
This sliding window is illustrated in Figure 3 by the 
coloured lines. This is equivalent to saying that disks 
are fading with time. Figure 6 shows size functions for 
the 0.4 < z < 0.6 and 0.6 < z < 0.8 bins with 0.2 and 0.4 
mag of fading. These models do little to the size func- 
tion at smaller sizes (/i < 4 kpc) but have some effect 
at larger sizes. About 0.2 mag is required to have the 
0.6 < z < 0.8 size function match the 0.8 < z < 1.0 
one, and 0.4 mag is needed from 0.4 < z < 0.6 to 
0.8 < z < 1.0. 

The 0.2 < z < 0.4 size function is shown in Figure 7 
with luminosity evolution ranging from 0.0 to 1.5 mag. 
Although neither the 0.5 mag nor no-evolution models 
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can match this size function to the one at high redshift, 
about 1.0 and 1.5 mag appear to be adequate depending 
on disk size. The 1.5 mag model is a better model for 
sizes < 4 kpc, and the 1.0 mag model is more reasonable 
at larger sizes. It is therefore possible that the largest 
galaxies may undergo a milder evolution than smaller 
galaxies. 

0.01 I 1 . 1 . 1 . 1 
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Fig. 6. — Size function of disk galaxies with luminosity evolution. 
The red dashed line shows the size function at 0.8 < z < 1.0 bin, 
and it is reproduced in all other redshift bins. The green and 
blue lines are the size functions at the indicated redshift with 0.2 
and 0.4 mag of luminosity evolution respectively. Their selection 
windows are shown in Figure 3. Black points are the no-cvolution 
size functions. 
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Fig. 7. — Size function of low-z disk galaxies with pure luminosity 
evolution. The coloured lines are the 0.2 < z < 0.4 size function 
with 0.5 (green), 1.0 (blue) and 1.5 (magenta) mag of evolution. 
Error bars for these models are the same as the no-evolution data 
points. The selection window for each size function is shown in 
Figure 3 with the same corresponding colours. Black points are the 
the no-evolution 0.2 < z < 0.4 size function, and the red triangles 
are the 0.8 < z < 1.0 no-evolution size function. 



6.3. Size Function from the Sloan Digital Sky Survey 

Studies of disks at high redshift can now greatly bene- 
fit from our increased knowledge of local disk properties 
thanks to the Sloan Digital Sky Survey (SDSS). A large 
dataset such as SDSS provides the ideal local baseline 
for comparison. We computed disk size functions from 
SDSS to compare with our CFHTLS results. The data 
come from Data Release Six (Adelman-McCarthy et al. 
2007) and the details of the SDSS size function analysis 
are given in Simard (2008, in preparation). Briefly, the 
three main galaxy selection criteria are Petrosian mag- 
nitude rp^corr (corrected for Galactic extinction), red- 
shift and spectral classification. We selected objects with 
14.0 < rp^corr < 17.7, 0.005 < z < 0.2 and the spectrum 
of a galaxy as defined by the keyword SpecClass in the 
SpecPhoto database table (SpecClass=2). These crite- 
ria yielded 522,453 galaxies. The nominal r— band sur- 
face brightness limit of the SDSS spectroscopic sample is 
/i50 = 24.5 mag arcsec"^ (Strauss et al. 2002). However, 
we set the faint surface brightness limit of our sample 
to /i5o = 23.0 mag arcsec"^ following Shen et al. (2003) 
to retain a complete sample. The final total number of 
objects satisfying our selection criteria was 493,366. The 
redshift distribution of this SDSS subsample peaks at 
z ~ 0.1. 

Galaxy structural parameters were measured from 
bulge-f disk decompositions performed using the G1M2D 
version 3.1 software package (Simard et al. 2002). We 
used the sum of an exponential disk and a de Vaucouleurs 
bulge (Sersic index n = 4) as our galaxy image model. 
Fits in (7, r and i were done using the separate fitting 
procedure described in Simard et al. (2002). GIM2D fit- 
ting failed for 226 objects (failure rate of 0.046%), but 
these "objects" were artifacts (e.g., image defects, bright 
star spikes) in the SDSS catalog. Galfit was not run on 
the SDSS images, but we were able to directly compare 
G1M2D/SDSS and Galfit/CFHTLS bulge fractions and 
disk scale lengths using 49 galaxies in the overlap region 
between the CFHTLS D3 field and the SDSS survey area. 
The comparison was done in the i-band because the typ- 
ical redshift of the SDSS galaxies is 0.1. The agreement 
is excellent as shown in Figure 8. 

The comparison between the SDSS and CFHTLS size 
functions is shown in Figure 9. The SDSS size function 
is based on the g- band fits for all redshifts which corre- 
sponds to the i-band at z = 0.9. Given that size changes 
as a function of colour are small (Harden et al. 2005; de 
Jong 1996), the effect of using a uniform filter for the re- 
spective CFHTLS and SDSS fits are not expected to have 
a significant effect on the size functions. The CFHTLS 
size functions are the no-evolution ones computed using 
the 0.8 < z < 1.0 luminosity-size selection window from 
Figure 3, and the no-evolution SDSS size function was 
computed using the same selection window. There is ex- 
cellent agreement between the 0.2 < z < 0.4 CFHTLS 
and SDSS {z ~ 0.1) size functions. This is a good check 
of the two independent size function calculations. The 
SDSS size functions with pure luminosity evolution were 
calculated by brightening the magnitudes of the SDSS 
galaxies before applying the CFHTLS selection window. 
This artificial brightening effectively selects the same ob- 
jects in the magnitude-size plane as the varying selection 
window described earlier. Galaxies with h < 4.5 kpc 
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appear to require about 1.0—1.5 mag of brightening be- 
tween z = 0.1 and z = 0.9 whereas 0.5 mag of brightening 
is a better match for the larger galaxies. This is generally 
the same trend seen in the CFHTLS data. The trend is 
more apparent in the SDSS size functions because they 
have much better statistics. The SDSS size function with 
pure number density evolution was calculated by simply 
multiplying the SDSS space densities by a factor of 3 be- 
fore applying the CFHTLS selection window. It provides 
a very good match to the 0.8 < 2: < 1.0 CFHTLS size 
function. 

7. DISCUSSION 

The fundamental variables of disk evolution are lu- 
minosity, size and number density. Some studies have 
measured changes in the zeropoint of the luminosity- 
size relation (i.e., surface brightness), and it is im- 
portant to emphasize that they could not disentangle 
size and luminosity evolution. Other studies provided 
stronger constraints by including number density infor- 
mation through the size function. The list of previous 
studies on disk evolution is given in Table 1. We divide 
the following discussion into pure luminosity and pure 
number density evolution for the sake of clarity with the 
understanding that in reality, both variables affect the 
population to some degree. 

7.1. Pure Luminosity Evolution 

The amount of luminosity evolution found by all pre- 
vious works ranges from to 1.6 mag. Simard et al. 
(1999) and Ravindranath et al. (2004) found little or no 
surface brightness evolution and pointed out that an ar- 
tificial evolutionary effect could be introduced by select- 
ing galaxies in different ranges of absolute magnitudes at 
different redshifts. This artificial effect was due to local 
galaxies not following a constant brightness relation (de 
Jong and Lacey 2000; Shen et al. 2003; Driver et al. 2005, 
Simard 2008, in prep.). Local disks with absolute mag- 
nitudes bright enough to bo visible at high redsliift are 
nearly 1 mag brighter in surface brightness than fainter 
disks. In section 6.1, we selected objects in the same 
manner as these groups and found that number density 
evolution was needed in order to reproduce observed size 
functions. We will revisit this in the next section. 

Lilly et al. (1998) and Sargent et al. (2007) found that 
the size function of disks was constant with redshift after 
1 mag of Imninosity evolution had been included. This is 
in agreement with what we found for the case of a varying 
selection window. Our own estimates are between 1.0 
and 1.5 mag with a weak dependence on size. As seen in 
Figure 7, the 1.5 mag evolution model seems to be too 
much for the largest disks in our sample {h > 4.5 kpc). 
The 1 mag model provides a better match for these larger 
disks. Our comparison between SDSS and CFHTLS size 
functions yielded a somewhat lower amount of evolution: 
1 mag for disks with h < 4.5 kpc and 0.5 mag for the 
larger ones. 

The luminosity evolution of a disk is obviously tied to 
the evolution of its constituent stellar populations. It is 
therefore interesting that even the largest amount of lu- 
minosity evolution is lower than expected from passive 
evolution of a stellar population. Passive evolution is 
thought to be responsible for at least 1.5 mag of evolu- 
tion since 2; = 1.3 in early type galaxies (van Dokkum 



and Stanford 2003), and perhaps more (Trcu ct al. 2005; 
Gebhardt et al. 2003). An estimate of 0.7 mag between 
0.0 < 2; < 0.7 from Cassata et al. (2007) is in rough 
agreement with our sample. However, the disks exam- 
ined here may have had a more active star formation in 
the past than the galaxies of these studies. This would 
suggest that our evolution estimate is likely lower than 
expected, although this amount is highly dependent on 
star formation history and the growth of stellar mass. 
The star formation history of the universe has declined 
sharply since z = 1 (Madau et al. 1996), and star for- 
mation rate is highly dependent on mass (Juneau et al. 
2005): less massive (and smaller) systems have been far 
more active in recent times than more massive ones. The 
pure luminosity evolution model here assumes that all 
disks have undergone the same amount of evolution. A 
greater luminosity evolution for smaller galaxies as seen 
in our data and previous works could be interpreted in 
terms of increased activity in lower mass systems (e.g., 
Zheng et al. 2007). 

7.2. Pure Number Density Evolution 

In section 6.1, we presented the size function of galaxies 
computed with a fixed selection window at all redshifts. 
The size functions at 0.2 < 2; < 0.4 and 0.8 < 2; < 1.0 
could be reconciled by increasing the number density of 
high redshift disks by a factor of 2.5 with no size evolu- 
tion (Figure 5) . A comparison with SDSS required an 
increase by a factor of 3 (Figure 9). Sargent et al. (2007) 
found that larger disks were less abundant at high red- 
shift albeit with a more modest increase in density; the 
largest disks in their sample {h >10 kpc) were only 60% 
as abundant as local disks. In determining what amount 
of number density evolution is actually needed for our 
sample, it is important to remember that we selected 
the disks of the most disk-dominated galaxies {B/T < 
0.2), and that disks that have recently grown a signifi- 
cant a bulge component would be excluded. The growth 
of 'classical' bulges is driven by mergers, whereas the 
growth of 'pseudobulges' is driven by internal (secular) 
processes within the galaxy. 

Mergers and the quenching of star formation are re- 
sponsible for transferring galaxies from the so-called 
"blue cloud" where most of the disk-dominated galax- 
ies in our sample would reside to the "red sequence" of 
quiescent galaxies (Bell et al. 2006; Faber et al. 2007). 
Although disks are relatively fragile, they are not en- 
tirely destroyed in merger events. Disks may re-form in 
merger remnants if the merging galaxies were sufficiently 
gas-rich (Springer and Hernquist 2005; Robertson et al. 
2006; Naab et al. 2006; Brook et al. 2007). However, 
all merger remnants would all have bulge components 
that would force them out of our sample. The epoch of 
merger-driven galaxy evolution is thought to have peaked 
early (2 ~ 2—3) in the history of the Universe (Kormendy 
and Kennicutt 2004). The merger rate in recent times, 
particularly since z = 1, is thought to have been quite low 
(Patton et al. 2000; Bundy et al. 2004; Lin et al. 2004; 
Blanton 2006). More specifically, Patton et al. (2000) 
find that only 6.6% of present day galaxies with rest- 
frame B-band absolute magnitudes — 18 < Mb < —21 
(similar to the range probed by our sample) have had 
a major merger since 2 = 1. Galaxies in the blue cloud 
are thought to follow a relatively quiet evolution with less 
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Fig. 8. — Direct comparison between Galfit/CFHTLS and GIM2D/SDSS structural parameters. Left-hand panel: i-band bulge fraction. 
Right-hand panel: i-band disk scale length in arcscconds. All SDSS data are from Data Release 6. 
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smaller than the size of the data points. Left-hand panel: Pure luminosity evolution. The SDSS size functions with 0.0, 0.5, 1.0 and 1.5 
magnitude of evolution include 41,056, 59,060, 74,747 and 86,832 disks respectively. Right-hand panel: Pure number density evolution. All 
SDSS data are from Data Release 6. 



TABLE 1 

Comparison of surface brightness evolution of different works 



Work 


Number of disks 


Redshift Range 


fi evolution 


Forbes et al., 1996 


17 


0.20 < z < 0.84 


0.6 - 0.85 


Roche et al., 1998 


270=^ 


0.20 < z < 0.90 


0.94 ± 0.23 


Lilly et al., 1998 


341 


0.10 < z < 1.00 


0.8 ± 0.3 


Schade et al., 1996a 


110 


0.50 < 2 < 1.10 


1.6 ± 0.1 


Simard et al., 1999 


136 


0.00 < 2 < 1.00 


minimal 


Ravindranath et al., 2004 


1508 


0.20 < 2 < 1.25 


< 0.4 


Harden et al., 2005 


5506 


0.00 < 2 < 1.10 


0.99"± 0.06'" 


Sargent et al., 2007 


12000 


0.00 < 2 < 1.00 


1.0 


This Work 


65000 


0.20 < 2 < 1.00 


1.0 - 1.5 



''to 2 = 3.43 
l^V band 
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than 10% experiencing a major merger which would then 
place them on the red sequence (Blanton 2006). There- 
fore, such a small fraction of mergers cannot yield sig- 
nificant number density evolution. Although the merger 
rates for more massive galaxies (M» > 5 x lO^^M©) may 
be slightly higher (Bell et al. 2006), it is unclear how this 
would relate to our disk sample given that it applies to 
a wide range of morphologies and mass-to-light ratios. 

As noted above, mergers are not be the only mech- 
anism that can transform pure disk galaxies into more 
bulge dominated ones. As described in Bower et al. 
(2006), field galaxies can undergo secular evolution pro- 
cesses which may produce a significant change in the 
number density of disks since z = 1. Disk instabilities 
and bars are examples of dynamical drivers which would 
cause significant morphological transformation without 
external influence. In these cases, gas within the disk is 
redistributed towards the centre leading to enhanced star 
formation and bulge growth (Kormendy and Kennicutt 
2004). The resulting "pseudobulges" have dynamically 
cold structures similar to disks, and their host galaxies 
arc typically located within the blue cloud (Drory and 
Fisher 2007). The timescale for the formation of a bar 
and then a pseudobulge is a few gigaycars (Kormendy 
and Kennicutt 2004; Drory and Fisher 2007), and the 
time elapsed since z = 1 would be amply sufficient to 
allow this process to operate on the disk population. 

8. CONCLUSIONS 

Using deep, wide-field imaging from the CFHT Legacy 
Survey covering two square degrees, we perform two- 
dimensional, bulge+disk decompositions on all galaxies 
brighter than i' < 24.5. Biases are carefully characterized 
through artificial galaxy simulations. We select the disks 
of the most disk-dominated galaxies with bulge fraction 
less than 20% (6000 galaxies), and we compute disk size 
functions over the redshift range 0.2 < z < 1.0 using 
different selection windows to illustrate the importance 
of selection effects. We also compute the size function of 
a very large sample of disks from the Sloan Digital Sky 
Survey to use as a local {z ~ 0.1) comparison. We find 
that: 

1 . The CFHTLS size functions computed by using the 
same fixed luminosity-size selection window at all red- 
shifts exhibit evolution that appears to be best modelled 
by a pure number density evolution. The z = 0.3 size 
function is an excellent match to the z — 0.9 one if disks 
at the highest redshift arc a factor of 2.5 more abundant 
than in the local universe. The SDSS size function would 
also match the z = 0.9 CFHTLS size function very well 
with a similar change in number density. However, cur- 
rent estimates of the galaxy merger rates below z = I 
are too low to produce such a significant change in the 
number density of the disk population in our sample. 

2. The CFHTLS size functions computed by using 
a varying luminosity-size selection window with redshift 
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remain constant if the selection window is shifted by 
1.0—1.5 mag towards fainter magnitudes with decreasing 
redshift. There is a weak dependence on disk scale length 
with smaller {h < i kpc) disks requiring more luminosity 
evolution than larger ones. This amount of pure lumi- 
nosity evolution is consistent with previous works and 
expectations from the evolution of disk constituent stel- 
lar populations. 

The evolution of disks is likely to be a more compli- 
cated interplay between size, luminosity and number den- 
sity than the simple pure luminosity and pure number 
density evolution models explored here. Models of disk 
formation and growth are still in their infancy. Obser- 
vational results are beginning to provide interesting con- 
straints, but more work remains to be done. Direct mea- 
surements of disk spectral energy distributions through 
full multi-color bulge+disk decompositions over a wide 
range of wavelengths would go a long way in furthering 
this effort as it would disentangle the evolution of disks as 
a function of morphology, mass-to-light ratio and stellar 
mass. 
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